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. Abstract 

, We study the non-equilibrium dynamics of a symmetry restoring phase transition 

on : in a scalar field theory, the "system" , linearly coupled to another scalar field taken as 

a "heat bath". The "system" is initially in an ordered low temperature phase, and 

the heat bath is at a temperature close to the critical temperature for the system. We 

CN . estimate the time at which the phase transition to the disordered (symmetric) phase 

C — ■ 

■ takes place. We derive, and integrate the one-loop effective equations of motion for the 
order parameter that include the effects of the heat bath. A semiclassical Langevin 

■ equation is derived and it is found that it contains a non-dissipative, non-Markovian 
kernel, the noise term is colored and correlated on time scales determined by the 
temperature of the heat bath. The range of validity of the Langevin equation and a 

O ■ consistent procedure to incorporate corrections are discussed. 

<D : 

& : 

; 1 Introduction and Motivation 

X- 

Non-equilibrium aspects of quantum field theory and statistical mechanics are beginning 
to receive considerable attention, for example within the context of structure formation in 
inflationary cosmology 0, 0, [| |j, dissipation and decoherence in quantum cosmology [|, in 
the theory of quantum brownian motion]?], |8|, |9|, and field theory in noisy environments [|I"0| 1 
to cite but a few. 

An important setting in which statistical mechanics out of equilibrium will play a very 
important role is in the description of the dynamics of phase transitions in the early universe. 
Recently the effective evolution equations to one loop order for the order parameter ITT] 



and the process of domain formation and growth Ji"2"| during a typical second order phase 
transition had been studied. 
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In this article we study the non-equilibrium aspects of a symmetry restoring phase transi- 
tion induced by coupling a system originally at low temperature in a broken symmetry state 
to a heat bath at a temperature close to the critical temperature of the system. Furthermore 
we derive a Langevin equation for the system that incorporates the effects of the heat bath, 
determine its range of validity, and provide a consistent procedure to incorporate corrections 
to the Langevin description. 

Phenomenological Langevin equations had been recently used in a variety of field theoret- 
ical applications: for a semiclassical non-equilibrium analysis of thermal activation fl3|, [IJJ, 
and the dynamics of kink-pair production in a heat bath||15||. These phenomenological 
Langevin equations incorporate a linear friction term and an uncorrelated (white) noise 
that satisfy the classical limit of the fluctuation-dissipation theorem. 

Our motivations for this article are twofold. First we want to study the non-equilibrium 
dynamics of a phase transition induced by coupling a system to a heat bath from a first- 
principle calculation. Secondly we want to qualify and quantify the situations in which a 
phenomenological Langevin equation with friction and white noise is justified, again, by 
establishing a first-principle derivation of the field-theoretical Langevin equation, its range 
of validity and possible corrections. One of the physically relevant questions that we want 
to address is about the time scales in which a "system" equilibrates with a heat-bath. 

An equilibrium description of a thermodynamic system assumes that the degrees of free- 
dom under consideration had reached thermal equilibrium with a "heat bath" or reservoir, 
and that the system is observed at times much larger than its typical relaxation time. When 
all the dynamical transient effects have had enough time to relax, the final state of the sys- 
tem would be that of equilibrium and the final state will not be sensitive to details of initial 
conditions or correlations, but only to the global thermodynamic properties of the "heat 
bath", like temperature etc. 

Although these assumptions are commonly used and clearly correct in a wide variety 
of experimental situations, their extrapolation to the realm of a Quantum Field Theory in 
extreme environments like the Early Universe or Heavy Ion Collisions is not only unclear 
but perhaps unwarranted a-priori. In particular, the concept of a "heat bath" is somewhat 
obscure when studying a field theory as a "system". Strictly speaking, a "heat bath" or 
reservoir, must have certain characteristics that allow to separate in a well defined manner 
the reservoir from the system, in particular the characteristic relaxation times of the heat 
bath must be much shorter than those of the system, secondly, the reservoir must have an 
infinite specific heat (a condition necessary for the equivalence between a micro canonical 
and canonical description) and thus infinitely many degrees of freedom. 
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The second condition is fulfilled if the heat bath is a quantum field theory, but if the 
"system" under consideration is also a quantum field theory, both system and bath have 
infinitely many degrees of freedom and the separation becomes somewhat artificial. 

The separation of time scales presents another conceptual difficulty in field theory because 
both system and bath have degrees of freedom at all energy scales with typically comparable 
density of states at large energies (determined solely by the relativistic dispersion relation). 

If the relaxation time for the heat bath is much shorter than that of the system, one may 
ignore the dynamics of the heat-bath and only ask questions about the system. Effectively 
this is achieved by "tracing over" the degrees of freedom of the bath. 

In this article we examine the situation in which a self interacting scalar field is coupled 
to another scalar field at a particular time. Prior to this time both fields are uncoupled and 
assumed to be described by a thermal equilibrium ensemble, but at different temperatures. 
This initial condition is somewhat arbitrary and unjustified, but our purpose is to study the 
time evolution of one of the fields (the self-interacting scalar field) taken to be the "system", 
and trace out the other field, taken to be the "heat bath". If the "system" field is at zero 
temperature in a state of broken symmetry, how long does it take to produce the phase 
transition and to restore the symmetry if the "heat bath" is at a temperature close to the 
critical??. 

A very succesful approach to study a Brownian particle in a heat bath or environment 
is provided by the Langevin equation. The effect of the heat bath is to introduce a friction 
term and a stochastic noise (typically assumed to be uncorrelated) both being related by the 
fluctuation-dissipation theorem. 

We want to study if such a description is suitable for a field theory and if so under which 
circumstances, establishing the range of validity of such description and possible corrections. 

In section II the essential elements to study field theory out of equilibrium are introduced. 
In section III we study the dynamics of the symmetry restoration phase transition induced 
by the heat bath. Section IV presents the derivation of the Langevin equation, discusses the 
range of validity and the formal procedure to incorporate corrections. Our conclusions are 
summarized in section V. Two appendices present many technical details. 

2 Evolution out of Equilibrium: Real Time Analysis 

As mentioned in the introduction, our aim is to study the time evolution of the "system" 
when it is coupled to a "heat bath" at a different temperature. By "system" we now refer to 
the simplified case of a self-interacting scalar field, and by "bath" we will take for simplicity 
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the case of a free massive field with a mass much larger than the mass of the quanta of the 
"system", so as to justify that the time scales of the heat bath are much shorter than those 
of the system. Eventually we will take the temperature of the heat bath to be much larger 
than the masses (of both system and bath). In this situation, we can relax the constraint 
that the mass of the particles in the heat bath be much larger than that of the system. 

One can avoid having to face the question of how the "bath" acquired its temperature, 
by assuming that the bath degrees of freedom are in a highly excited state corresponding to 
a microcanonical ensemble with an energy density £ ^> m a with m a the mass of the particles 
in the bath. Under these circumstances one can pass from the microcanonical ensemble at 
this energy density to a canonical ensemble by a Laplace transform. The temperature of 
the equivalent canonical ensemble is T a w £*. We will use the description in terms of the 
canonical ensemble. 

We will also require that the coupling of the system to the bath satisfies the following 
conditions: i) locality ii) manifest Lorentz covariance, iii) maintains the renormalizability of 
the scalar field theory. 

For calculational simplicity we will also assume that the coupling is linear in both fields. 
The simplification of a "free field bath" and linear coupling will allow us to present detailed 
calculations, these restrictions may be relaxed at the expense of further complications. 

Thus our situation may be modelled by introducing the time dependent Hamiltonian 

H(t) = H$ + H a + Hj(t) (1) 

where 

H* = I {^niOr) + \{V<$>{x)f + \(-n 2 )$ 2 (x) + ^ 4 (x)| (2) 

is the Hamiltonian of the self-interacting field that will be taken as the "system" , and /i 2 is 
chosen to be positive with the purpose of studying non-equilibrium aspects of phase transi- 
tions. 

H„ = J v d s x[lnl(x) + \{Va{x)Y + \mla\x)) (3) 

is the Hamiltonian for the free field taken to be as the "bath" . Finally Hj(t) is the interaction 
Hamiltonian given by 

Hj(t) = gG(t) f d 3 x<S>(x)a(x) (4) 
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To completely determine the time evolution of the "system- bath", we need to specify 
the initial conditions. We assume (without any a-priori justification), that for t < both 
the system and bath are in thermal equilibrium at temperatures T$ ; T CT respectively. Fur- 
thermore we will assume that the field $ is in a broken symmetry state at T$ « and 
that T a ^> m^,/i 2 . Notice that the coupling constant g has dimensions of (mass) 2 , and a 
perturbative expansion in this coupling, will involve the ratio of this coupling to a particu- 
lar energy scale of the heat bath. Under the assumption that the temperature of the bath 
is much larger than the masses of the bath and system field, we expect the perturbative 
expansion parameter to actually be in terms of the dimensionless ratio g = g/T 2 . Further- 
more, in this limit, we expect the characteristic time (relaxation time) of the heat-bath to 
be t bath « 1/T a . 

For the study of the dynamics in real time, the quantity of interest is the density matrix 
in the Schroedinger picture. For t < the system and bath are uncoupled and the initial 
density matrix is given by 

p(0) = e ~P* H * ® e~ p " H " (5) 

where (3=1/ KbT for the respective temperatures. In the Schroedinger picture, the density 
operator evolves in time as 



p(t) = U{t)p{Q)U-\t) (6) 

with U (t) the time evolution operator. 

As mentioned previously, our aim is to understand the dynamics of a field interacting 
with a heat bath. Physically we are interested in correlation functions of the scalar field $, 
but not on any properties of the heat bath. Thus, we proceed to trace over the degrees of 
freedom of the heat bath, obtaining an effective theory for the scalar field. 

This procedure, starting from a closed "system-heat bath" and tracing over the reservoir 
degrees of freedom leads to the description of the "system" as an open environment, and is 
fundamentally equivalent to passing from a microcanonical to a canonical description. 

The procedure of tracing out (performing the path integral) the bath degrees of freedom 
yields to a "reduced density matrix" for the system. Although the full density matrix 
satisfies the Liouville equation with the full Hamiltonian, the reduced density matrix satisfies 
a Fokker-Planck equation [0, pTTf] . 

This procedure originally proposed by Feynman and Vernon |18[|, has been used re- 
peateadly in the literature |TB|, in particular in the study of dissipative effects in quantum 



mechanics and for the consistent treatment of quantum Brownian motion fl^. There are 
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some major differences between our work and past work on quantum brownian motion and 
dissipative systems. In the usual treatment of quantum brownian motion [0, a particular 
density of states must be assumed for the heat bath. The choice of density of states leads to 
different dissipative behavior, and is reflected in the correlation functions of the heat bath 
through the fluctuation-dissipation relation. 

In our case, the nature of the heat-bath, and its coupling to the system, are constrained 
by the requirements of locality, relativistic covariance, and renormalizability, which severely 
constrain the density of states of the heat bath and the possible couplings to the system. 

The expectation value of an arbitrary operator (0)(t) is given by 

It is more illuminating and convenient for our purposes to cast the above expression in 
terms of the time evolution operator U. This is achieved by first choosing a large negative 
time T < (not to be confused with temperature), for which U(T) = U$(T)U a (T) where 
Ui(T) = exp[— iH(T\ (for i = $ ; a) we may write exp[— (3iHj\ = exp[—iHi(T — ifa — T)] = 
Ui(T-iPi,T). 

Inserting in the trace U (T)U(T) = 1, commuting f/ _1 (T) with p(0) and using the 
composition property of the evolution operator, the expectation value (|7|) becomes 

= Tr {U*(T - ijh, T)U a (T - i(3 ai T)U~\T, t)OU(t, T)} 
WU Tr {U$(T — i(3$, T)U a (T — T)} U 

It proves convenient to choose a large positive time T' and to insert U(t,T')U(T',t) = 1 
(with U the full time evolution operator) to the left of O in (|]) to extend the contour in the 
numerator to an arbitrary large positive time T"|[20||. Finally, the thermal expectation value 
of the operator O is 

= Tr{{U*(T - ijh, T)U a (T - ifc, T)}U(T, T')U(T', t)OU(t, T)} 
[ ){) Tr{U 9 (T-il3 9i T)U ff {T-il3 a ,T)} U 

In the complex time plane, the numerator represents the following process: from T < 0, 
evolve in time up to t, insert the operator O, evolve further up to T', and backwards from 
T' to T, finally down the negative imaginary axis to T — i(3$ for $ field (T — if3 a for the a 
field) depicted in figure (1). The denominator just evolves along the negative axis from T to 
T — 2/3$ ;(7 for $ ; a. Eventually we will take the arbitrary times T — ► — oo ; T' — >• oo. 

The insertion of the operator O may be obtained as usual by introducing external sources 
coupled to the particular operator. 
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Since we propose to use perturbation theory to study the non-equilibrium correlation 
functions, it proves more convenient to introduce different source terms for all the time 
evolution operators in the trace. The sources will be different on the different branches 
along the complex time contour. However, we do not introduce source terms for the field a 
as it will be traced out, and no correlation functions of this field will be computed. We are 
thus led to consider the following generating functional 



Z[J£, Ji, 4} = Tr{[U*(T - i/3», T; jg)U a {T - T)}U{T, T'; J^)U(T', T; J+)} (f 0) 

where J$ are the source terms for the field $ along the forward (+) and backward (— ) 
segments of the contour, and the denominator in (H) is given by Z[0, 0, J 13 ]. This formalism 
allows to study more general situations with different choices of p(0). 

The necessity for these generating functionals to study non-equilibrium quantum statis- 
tical mechanics was originally proposed by Keldysh|22| and Schwinger ||2f|| and has been used 
in the literature very 



often |2|, H, H, H, |7|, ||, The trace over the a field (functional integral) may be 



easily performed because it is a free field linearly coupled to the scalar field $, which acts as 
a "source term" for the a field, yielding to a non-local quadratic influence functional§\T 



- l -g 2 J o dt x J o dt 2 $\t x )Df{t 1: t 2 )$ a {t 2 )\ ill) 

with a sum over the indices a,b — +, — , corresponding to the forward T — > X" and backward 
T' — >• T branches of the contour respectively, and (D% b (ti, t 2 )) are the Green's functions on 
the contour for the "heat-bath" given in appendix A. 

We relegate most of the technical details to the appendices, where we show that the final 
form of the generating functional becomes 



r j roo poo ~i 

exp / dtj dhJ^DtiUM)^)) (12) 

with a, b = +, — . 

The Green's functions on the contour and their properties are analyzed in appendix A. 
We are now in condition to study the dynamics in real time for the field $ in the presence 
of the "heat-bath" . As mentioned in the introduction, one of our principal motivations is to 
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address the fundamental question of symmetry restoration by the heat bath. In particular, if 
the system is originally at very low temperatures (or zero temperature) in a broken symmetry 
state and is suddenly coupled (at t — 0) to a "heat bath" in equilibrium at a temperature 
close to the critical temperature, how long does it take for the symmetry to be restored?? 
What are the relevant time scales for the dynamics of symmetry restoration?? The answer to 
the last question is certainly far from obvious. There are several widely different time scales 
in the problem, the typical time of the heat bath (determined by the mass and temperature 
of the a field), the mass and temperature of the scalar field (system), and finally the coupling 
system-bath. It is not a-priori obvious which time scale or combination thereof will determine 
the dynamics of the process of symmetry restoration. However, on physical grounds, we 
expect that for very weak coupling between the system and the heat bath, this time will be 
very large. This intuitive argument must be qualified however. The dimensionless coupling 
is expected to be g/T%, and the transition time is expected to be "large" on time scales of 
the heat bath. Thus when the temperature of the bath is much larger than the masses of 
both system and bath we expect the "transition time" t c (the time at which the symmetry 
restoration phase transition occurs) to be such that 

t c T a > 1 

3 Symmetry Restoration 

Usually the physics of a phase transition is studied by means of the static effective potential 
(free energy density), as the position of the minima of the effective potential determine the 
thermodynamic equilibrium states of the system. If the minima correspond to a nonzero 
value of the expectation value (or thermal average) of the scalar field the symmetry is 
spontaneously broken, otherwise the symmetry is restored. The vacuum expectation value 
or thermal ensemble of the scalar field serves as the order parameter for the transition. In a 
translational invariant theory this order parameter is independent of the spatial coordinates, 
and in equilibrium it must be independent of time. 

In our case, we couple suddenly the system to the heat-bath at time t — and let the 
system evolve in time with the full interaction Hamiltonian. Clearly this is not an equilibrium 
situation, as the initial density matrix does not commute with the full Hamiltonian for t > 0. 
The static effective potential is not the proper quantity to study the dynamics of the phase 
transition and one should invoke the effective action to account for the real time dependence 
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of the situation. Because translational invariance is still preserved, the order parameter 

only depends on time (and certainly temperature). 

Initially the scalar field is in equilibrium in the broken symmetry phase at very low tem- 



perature (T$ <C T c = y24/i 2 / A) and is suddenly coupled to the heat bath whose temperature 
is T a f» T c 3> m<j , \i . The order parameter will evolve in time and the phase transition will 
occur when the order parameter becomes zero or begins to oscillate around zero. 

Obtaining the effective action even to one-loop order is clearly an imposing task, the same 
information about the phase transitions is obtained from the effective evolution equations for 
the order parameter. These equations of motion are formally obtained as a variational princi- 
ple from the effective action, by requiring that the functional derivative of the effective action 
with respect to the order parameter vanishes. It turns out that it is much easier to obtain the 
effective equations of motion for the order parameter from the non-equilibrium formalism. 
This is achieved by using the tadpole method and applied to this non-equilibrium situ- 
ation by shifting the fields $ ± on the forward and backward branches as followsfll], [T^lfsee 
appendix B) 

$±(r,t) = <f)(t) + ip ± (f,t) (14) 

The reason for shifting both (±) fields by the same configuration, is that enters in the time 
evolution operator as a background c-number field, and evolution forward and backwards 
are now considered in this background. The effective equations of motion are obtained by 
requiring that 

(V ± (r,t))=0 (15) 

These two conditions consistently result in the same evolution equation for the order param- 
eter as a consequence of the cyclic property of the trace. 

As explained in appendix B the evolution equations are obtained by expanding the action 
in ([51]) around <p(t), the quadratic terms for ip will define the propagators and the higher 
order terms will be treated in perturbation theory. To one loop order we find 

|^(t) + + ^ 3 (t) - f Q dt'sm[m a {t-t')W) 



+ 2^ ) /(2^ ( ^ (t) ^ (t)) = ° (16) 
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The retarded kernel in the above expression is arising from the influence functional, notice 
that it does not depend on the temperature of the heat bath, this is a consequence of the 
fluctuation- dissipation theorem. This kernel is non-Markovian (has memory), but perhaps 
more importantly it is non-dissipative. 

We will see later that this is the retarded kernel entering in the Langevin equation that 
describes the effective semiclassical evolution equations in the presence of the heat bath. 

The one-loop term (<^i(t)(/9 + -(t)) is recognized as the spatial Fourier transform of the 
equal time two-point correlation function of the fluctuations. This correlation function is 
obtained by inverting the operator of quadratic fluctuations. 

The quadratic terms that define these propagators, however, depend on the time de- 
pendent order parameter 4>(t), and also receive a contribution from the influence functional 
(that mixes ip^). The identification of the "free propagator" is extremely difficult as one 
must find the inverse of the time dependent quadratic operator for the fluctuations. This 
involves summing the Dyson series for the mixing term in the influence functional. Because 
of the complicated time dependence we were unable to sum up the Dyson series and have to 
content ourselves with a perturbative expansion for both the time dependence of the order 
parameter and the contribution of the influence functional. Thus we write 

-t* + ^(t) = M 2 + ^(0 2 (t) - 2 (O)) 

M 2 = V + ^ 2 (0) (17) 

with <p 2 (0) = 6/i 2 /A being the minimum of the tree level potential. We identify the first 
term on the right hand side of ( |TTD as the mass of quanta in the broken symmetry state. 
The second term will necessarily be of order \g 2 or higher, thus O(X) smaller than the term 
arising from the influence functional because the time dependence is induced by the coupling 
to the heat bath. As may be seen from (fl6|) , such a term will contribute to 0(X 2 g 2 ) to the 
evolution equation. We will consider only the first order correction in g 2 to the propagator 
arising from the influence functional, a typical diagram is depicted in figure 2(b). The one 
loop contributions to this order to (|T6| ) are depicted in Figure 2(c). Thus we will consistently 
neglect to this order the second term on the right hand side of fll7|) (time dependence of the 
order parameter). 

The first one- loop contribution depicted in Figure 2(c) is the familiar term, it is time 
independent and gives the usual result in terms of the temperature of the field for our 
purposes assumed to be 7$ rs 0. We will absorb this term in a renormalization of the bare 
parameters, fi 2 ; A. 
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Then we define consistently to this order the time dependent mass 



m 2 (t) = (-/4) + Am(t) (18) 

where the 0(g 2 ) contribution to the one-loop correlation function is depicted in figure 2(b). 
A simple analysis reveals, that for t ^ 0, this one loop diagram is ultraviolet finite, this 
contribution vanishes at t = 0. Thus there are no further renormalizations (up to this order) 
induced by the coupling to the heat bath. 
In the limit of T$ »s we find 



Am(t) = Am (0) (t) + Am (Tl) (i) (20) 



where 



A m (°)(t) = -ff 2 / , 2 _ / dt x / dt 2 sin[4(t-ti)] [sin[w fc (ti - t 2 )] cosK(t - t 2 )] 

Z ./ (271") UJfcUJfc Jo Jo 

+ cos[a> fc (ti-t 2 )]sin[a;j;(t-t 2 )]] (21) 

^k Tff) W = A^ 2 / /V f rft 2 {sinK(t - oosfofa - t 2 )] 

J (27r) d w^Jfc •'o jo 

x sinK(t-t 2 )]} (22) 



^ = \/k 2 + M 2 ; u k = \Jk* + ml (23) 

Let us first consider the finite temperature contribution. For this purpose it proves 
convenient to introduce the following dimensionless variables, 

x = ^r , r = T a t (24) 

After performing the time integrals in the high temperature limit (T a ^> m a ; M 2 ), the 
finite temperature contribution (p2|) becomes 



Am^(r) = ^ff^l fds , 1 J ^[l-cos(2xr)] + ^[l- Sin(2 ^ ]| (25) 
V; 2n 2 \T*)Jo x(e x - 1) \ 16x 2 L V ;J 8 1 ir J 

This expression clearly shows that in this limit, the time scale for the heat bath is 
1/T a . As argued previously, we expect that for small system-bath coupling (g/T 2 -C 1) the 
transition time t c ^> 1/T a . 
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For t 3> 1, the integrand in ( |25| ) is sharply peaked at x ~ allowing a saddle-point 
approximation to the integral. We find for r ^ 1 



Am^'fr) » 0.012r 3 AT 2 




(26) 



In the same limit we estimate the contribution coming from the zero-temperature part 
to be 



Am(°'(r) 



0.013r 2 AT 2 |j 



(27) 



Thus as expected, in the r>l limit, the finite temperature contribution dominates. 

It is straightforward to see that in this limit the system-bath dimensionless coupling is 
g 2 /T£, and that the terms that we neglected in a perturbative expansion are of the order 
0(X 2 , \ 2 (g 2 /T£) , ■ • ■). Recalling that the critical temperature for the scalar field $ is 
given by T 2 = 24/4/ A#, it is convenient to rescale the field and time to cast the equation in 
terms of dimensionless variables. Define 



s = fiRt ; x 



6/4 



(2? 



the equation of motion (|I(J) becomes 



d 2 X (s) 
ds 2 



+X 3 (s) 



X{s) 
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1 -0.29s 3 f 



m a fi R Jo 



ds' sin 



s — s 



X(s 



0.31s 2 ^ 







(29) 



Admittedly inserting the result from the one-loop term in the equation of motion is not 
quite consistent at early times because the one-loop results ( p6| , p7|) were obtained for large 
r = s(T a /fi[{). However, by integrating the above equation we will obtain a qualitative 
understanding of the time evolution of the order parameter. 

Roughly speaking, the phase transition from the broken symmetry state to the disordered 
phase takes place when the effective time dependent mass ( |T8D becomes zero. This would 
certainly be the case in the absence of the non-Markovian term (retarded kernel in flTB])). The 
effect of this term may be understood by expanding the field around the broken symmetry 
values 4>± = ±\j6[i R /\ji, and considering the deviation from these values as perturbations. 
One finds that the retarded kernel plays the role of a forcing term in the equation of motion 
for the deviation, that tends to shift the value of the minima further away from the origin. 

However, we still expect that the phase transition will occur at time scales when the 
effective mass term vanishes, as the oscillations of the field will now occur near zero. 
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From the above effective equation of motion we find that the effective mass ([18]) vanishes 
at the "critical time" t c given by 

i 

(30) 

This is one of the main results of this work. Although the time scale is determined by 
the temperature of the heat bath, the critical time is a non-analytic function of the coupling 
between the system and bath. 

Figure(3) shows the evolution of the equation of motion ( p9l) for the (arbitrary) value of 
the parameters llr = m a ; g = m 2 ; T a = T c = 10m a . For these values, the critical time 
predicted by (|30|) is fiRt c ~ 5.8, quite consistent with the time at which the field crosses zero 
in Figure (3). 

Thus we see that for T a « T c and g 2 /T* C 1, i c > 1 /T a (or alternatively t c 1) quite 
consistent with our original assumptions. 

Notice that the critical time is non-perturbative in terms of the system-bath coupling, in 
a sense very much like the critical temperature is non-perturbative in terms of the scalar field 
self-coupling. One would then argue that the higher order corrections will only contribute 
perturbatively to the critical time. This argument however ignores the possible infrared 
divergences arising near the phase transition as it happens in equilibrium finite temperature 
field theory near the critical point. So we expect that, very much as it happens in equilibrium 
finite temperature, approaching the critical point or "critical time" in a reliable manner will 
require a non-perturbative resummation of the infrared sensitive diagrams [[!(], |31], |32|. Thus 
just as the critical temperature signals the breakdown of perturbation theory to study the 
phase transition, this critical time may also signal the breakdown of perturbation theory 
in the system-bath coupling. Presumably a non-perturbative resummation technique would 
have to be invoked to study reliably the system-bath dynamics for times close to the critical 
time t c . 

This possibility would have to be studied further and is beyond the scope of this article. 

4 Semiclassical Langevin Equation 

A phenomenological but rather successful approach to study the non-equilibrium dynamics 
of a particle coupled to a heat bath in classical statistical mechanics is provided by the 
Langevin equation. This approach provides a satisfactory description of Brownian motion of 
a particle. The Langevin equation is the classical equation of motion modified phenomeno- 
logically by (basically two) terms that account for the interaction with the heat bath. A term 



T A \ T 2 

3 ' 5 # 5f 
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proportional to the velocity of the particle that incorporates "friction" and dissipation, and 
a stochastic force term that reflects the random "kicks" of the heat bath upon the Brownian 
particle. In most applications, this stochastic noise is assumed to be "white", that is com- 
pletelety uncorrelated. The coefficient in the friction term determines the relaxation time of 
the Brownian particle. The friction term and the stochastic force are ultimately responsible 
for the approach to equilibrium of the particle with the heat bath, and are thus related by 
the fluctuation-dissipation theorem. 

A very clear microscopic description leading to Langevin dynamics within the context of 



one-particle quantum mechanics has been presented by Caldeira and Leggett and Schmid||16 
p3| . Their study reveals that in general the friction term arises as a local approximation to a 
non-Markovian kernel for a particular choice for the density of states of the heat-bath, and 
that the noise term is typically correlated over time scales determined by the typical scales 
of the bath. 

An attempt to obtain a microscopic description of Langevin dynamics in field theory has 
been reported by Morikawa [j34j1 who considered a "fermionic bath". 



In this section we offer a semiclassical derivation of the Langevin equation for the case 
of a bosonic heat-bath, emphasizing its range of validity, and the formal steps to go beyond 
the semiclassical Langevin equation. 

Phenomenological Langevin equations have been recently used to study the dynamics of 



semiclassical configurations in a presence of a heat bath|T3|, [14], [L5[]. Typically these Langevin 
descriptions assume a friction term and a "white-noise" . We will show in what follows, that 
these assumptions may not be justified in many cases, and that a physically correct Langevin 
description must necessarily incorporate details of the particular "heat-bath" . 

We now consider the case in which the $ field is initially at zero temperature, and it is 
coupled at t = to the heat bath, again modelled by the free field a. The initial density 
matrix is now 

P(0) = (|0>»(0|#)®e-^ H - (31) 

Evolving this density matrix in time as in the previous section (details are provided in 
the appendices), and tracing over the bath variables (thus obtaining the influence functional) 
we obtain the following generating functional 

Z[Q] = J D$! J ^e l f-™ dixU[ * +] ~ u[ ^ ] } x (32) 

with the boundary conditions $ + (r, t = oo) = $~(r, t = oo) = 



14 



At this stage it proves convenient to introduce the "center-of-mass" and "relative" field 
coordinates \I/ and rj,^/ = i[$ + + <3> - ]; 77 = [$ + — $ - ]. These are recognized as the coordinates 
used in the Wigner transform of a density matrix [IS]. The boundary conditions on the fields 
now become \&(f, t = 00) = $ x ; r)(r, t = 00) = 0, the integral over $1, just defines the path 
integral over all field configurations, and finally the partition function (E32|) becomes 



Z[0] 



'•S[v,v] 



/OO ^ /-OO 

^ rf 4 x 2 6(ti - t 2 )/C/(ri, ti; f 2 , t 2 )tf (r 2 , i 2 ) - - y_ d 4 zi*(xi)^ 3 (xa) + 



0* /"CO /*oc 

z— / d 4 Xx / d i x 2 r](x 1 )IC R (f l ,t 1 ;f2,t2)r](x2) 
Z J— 00 J —00 



(33) 



where the kernels /Cr, /Cj are spatially translationally invariant and their spatial Fourier 
transforms are given by (see appendix A). 



K.i(k,t 1 ,t 2 ) 



— cos[cD fc (ti - t 2 )][l + 2n a (u k )]G(ti)e(t 2 ) 

ZU k 

—sm[uj k {t 1 -t2)]e{t 1 )e{t 2 ) 



2u k 



(34) 



with Cbk being the frequencies of the heat-bath. These kernels arise from the two-point 
correlation function of the heat-bath fields, the real Kr and imaginary part K,j are related 
by the fluctuation-dissipation theorem. 

To make contact with Langevin dynamics, it becomes convenient to cast the quadratic 
term for the r\ field as resulting from a Gaussian integral over a stochastic noise term with 
probability distribution given by 



V[£] = exp 

Now the final form for the partition function becomes 

Z[0]= \vi J V^Vr)V[£\e iS *ff [ *^} 
with the effective action S e ff given by 



(35) 



(36) 



Seff[*,V,€\ 



d A X\iq{x\ 



:-U-{-y?))*{ Xl )-^\ Xl ) 
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-2g 2 J d*x 2 Q(t x - t 2 )Kj[f u h; r 2 , t 2 )*{x 2 ) + Z{ Xl ) 

- d 4 x^{x)r 1 \x) (37) 

A semiclassical approximation to the above partition function requires the configurations 
that extremize the effective action. In particular the condition SS e ff/Srj = leads to the 
lowest order semiclassical result 

(□ + (-/i 2 ))*(n,*i) + A^3 (fi;ti) + 

2g 2 J dt 2 d 3 r 2 e(tt -t 2 )/C / (f 1 ,t 1 ;f 2 ,t 2 )^(f 2 ,t 2 ) = £(^0(38) 
with the Gaussian noise correlation function 

£i)£(r 2) t 2 )} = sfK R (r u h; r 2 , t 2 ) (39) 

This is the typical Langevin equation. This equation, however, is not the only dynamical 
evolution equation in the semiclassical limit, other equations are obtained by performing the 
variational derivatives with respect to the "noise" term and Langevin dynamics results 
from the decoupling approximation to these equations, in particular neglecting the coupling 
between \1/ and r\ and r\ and the "noise". A consistent improvement over the semiclassical 
Langevin equation will involve a perturbative expansion in these terms. Allowing for the 
couplings of the different fields to the "noise" term will introduce corrections to the "noise" 
correlation functions. This corrections are a manifestation of the "back-reaction" of the 
system on the bath correlations. 

One may obtain an equation of motion similar to the effective equation found in the 
previous section for the order parameter (|Tj]) from the above Langevin equation by splitting 
the field as 

*(f,t)=<Kt) + <p(f,t) (40) 

requiring that the "mean field" <p(t) obeys a source free equation and considering (p as the 
fluctuation whose equation of motion contains the noise term. The linearized equation for 
the fluctuation with the noise term, and non-Markovian kernel may be solved by introducing 
the retarded propagator. In the equation for the "mean-field" there is a term of the form 

^(t)(^(ri,t)) 2 
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By taking the average of this term over the noise with the noise auto-correlation function 
(|39l) one obtains the 0(\g 2 ) correction to the equation of motion of the order parameter, but 
not the one-loop correction without the bath (0(A)). This semiclassical Langevin equation 
does not incorporate quantum loop corrections of the scalar self-interaction. 

There are two features of the semiclassical Langevin equation that deserve comment: the 
first is that the non-Markovian kernel does not lead to dissipation. Secondly, the noise does 
not have "white" (delta functions) correlations. In fact these two features are related by the 
fluctuation-dissipation theorem. 

The fact that there is no dissipation is a consequence of the simplicity of the model, 
one could try a coupling to the heat-bath that is linear in the scalar field and quadratic 
in terms of the bath fields. Integration over the bath fields will also yield to a non-local 
influence functional which however will not be quadratic in the scalar field, but may be 
studied perturbatively. The quadratic contribution of the scalar field will involve a "bubble" 
diagram from the bath fields. This diagram will have a two-particle threshold that may 
lead to dissipative processes, which however, is very high in energy if the bath particles are 
much heavier than the system particles. Thus in this case long wavelength low frequency 
components of the field will evolve without dissipation. 

Such a diagram appears naturally in the case in which the bath field corresponds to 
fermions coupled to the scalar field via a Yukawa coupling. A memoriless (Markovian) 
friction term in the Langevin equation (proportional to the time derivative of the field) may 
be obtained at low frequencies provided the non-local kernel has a power series expansion 
in the transferred frequency with a linear term in the frequency. The coefficient of this 
linear term will be the friction coefficient. From the above discussion, it is clear that this 
possibility is not very easy to achieve if the bath particles are more massive than the system 
particles. In this case the multiparticle threshold will be very high in frequency and there 
will be no dissipation for low transferred frequency. This may be inferred from the work 
of Morikawa||34||. This situation may change if there are collective excitations in the bath, 
that produce a linear frequency dependence of the non-local kernel at low frequencies. This 
scenario, however, is not generic and will depend on the details of the heat bath. 

This conclusion is particularly meaningful within the context of Langevin dynamics of 
semiclassical field configurations. These are large amplitude coherent configurations, but 
typically slowly varying in space and time, thus corresponding to small frequency and mo- 
mentum transfers to the bath. Although very massive, because they are mainly composed 
of long-wavelength, low frequency modes, it is difficult for these configurations to "decay" 
and thus dissipate. 
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Derivative couplings between system and bath fields may lead to dissipative terms, but 
are very dangerous from the point of view of renormalizability. Couplings that are quadratic 
(or higher order polynomials) in the scalar field will yield to higher order polynomials in the 
influence functional and certainly will not lead to a simple Langevin description. 

Thus in order to obtain a semiclassical Langevin description there is a strong restriction 
on the system-bath coupling: the system must couple linearly to the bath, so that the 
influence functional is quadratic. A friction term may arise if the bath degrees of freedom 



are lighter than the system's degrees of freedom, but as discussed by Morikawa||34|| this will 
have to be studied case by case, in particular it may require a higher order calculation. 

The point of this discussion is that writing down a simple Langevin equation with a local 
dissipative term and a noise term with white correlations may have to be justified by looking 
at the particular models in detail. 



5 Conclusions: 

We have studied the dynamics of a phase transition in which a scalar field theory (the 
system) , initially in the ordered (broken symmetry) phase at low temperatures is coupled 
to a heat bath. This heat-bath is represented by another scalar field, linearly coupled to the 
"system field" and whose temperature is close to the critical temperature for the system. We 
derived, the one-loop effective equations of motion for the order parameter that incorporates 
the effects of the heat bath. 

The heat bath introduces a non-Markovian, non-dissipative kernel that is temperature 
independent and also introduces one-loop corrections that are temperature dependent. 

The effective evolution equations were integrated and an estimate of the time to complete 
the phase transition to the symmetric phase was obtained. This time is a non-analytic 
function of the "system-bath" coupling and may signal the breakdown of perturbation theory 
in terms of this coupling. 

We derived a semiclassical Langevin equation for the system, it contains a non-Markovian, 
non-dissipative kernel, and the noise term, although Gaussian (in this approximation), is 
colored and correlated on time scales of the order of the inverse temperature of the heat-bath, 
when this temperature is much larger than the masses of the fields. It is pointed out that the 
non-Markovian kernel will be model dependent and may only be approximated by a local 
dissipative (friction) term only in very special cases. The derivation of the Langevin equation 
permits to identify a formal expansion to improve over the semiclassical approximation and 
to account for quantum effects (loops) of the system field, as well as for back reaction of the 
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system on the bath. 
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Appendix A : Finite Temperature Generating Functionals 

In this appendix, we summarize the most relevant technical details leading to the gen- 
erating functional flT2"|) in section 2. Most of the steps may be found in the literature 



26| , [27], go[ |29| , however we present here the generalizations and modifications appropriate 
to our case in order to make our article self-contained. 

Starting from the generating functional ( |lOl) let us insert the resolution of the identity in 
terms of a complete set of field eigenstates 

J D<$>Da\<$>,a)(a,<$>\ = 1 (41) 

between all the time evolution operators, obtaining 

Z[J£, J$, 4} = J D$ 1J D$ 2 D$ 3 J DaxDa 2 Da 3 

{($!, (n\U*(T - z/3*, T; 4)U a (T - ij3 a , T)|$ 2 , a 2 )} x 
{($ 2 ,a 2 |f/(T,T / ;J^)|$3^3)($3,a 3 |f/(T / ,T; 4)1$!,^)} 

(42) 

Each matrix element in the above expression has a functional integral representation. Then, 
the generating function becomes (to avoid cluttering of notation we suppress the space-time 
indices) 

Z[J£,JZ,J%] = J D^D^D^ J V§ + V§-V^ 

| e i/ T T 'r*[*+]+4* + -^[*-]-^*-| je'/r^^l+^J x ^[$+,$-] 

(43) 

with the boundary conditions $+(T) = $^(T - = $i;$~(T) = $^(T) = $ 2 and 
$ + (T') = $~(T") = $3. This can be regarded as a path integral along the contour (C) in 
the complex time plane shown in Figure (1) with periodic boundary conditions. 

The non-local functional T is the "influence functional" WS, [T?J obtained by tracing-out 
the bath degrees of freedom [a) 

^[$+,$-] = J Da 1 Da 2 Da z J Va + Va-Va^ 
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with the boundary conditions u + (T) = a^{T— i/3 a ) = $1; cr~(T) = cr^(T) = a 2 and a + {T') = 
a~(T') = a 3 . The normalization term JF[0, 0] will cancel against the denominator in the 
computation of correlation functions. The influence functional can also be recognized as a 
path integral along the contour (C) shown in Figure (1). 

For the linear coupling between system and bath fields, the trace (path integral) over the 
bath variables can be done at once leading to the expression (|TT1). 

The propagators D^ b (a, b = +, — ) are constructed from the functions D>, D< which are 
the homogeneous solutions of the quadratic form in L a obeying the Kubo-Martin-Schwinger 
boundary condition (periodicity in imaginary time) 

D>{r 1 ,t l -i^-r 2) t 2 ) = D<{f l) t 1 \r 2 M) (45) 
The spatial Fourier transforms of these functions are given by 



iD>(k;t,t') = -L{(l-n CT K))e-^( t - t ')+n a (d; fc )e^^')} 

D>(k;t,t') = D<(k,t',t) 
1 



n °( Q k) = _ 1 ; oj h = \Jk 2 + ml (46) 

The Green's functions (D® b (ti,t 2 )) that enter in the integral are now explicitly given by 



D+ + (fi, tr, f 2 , t 2 ) = D>(? u t x \ f 2 , t a )6(*i - t 2 ) + D<(n, h; f 2 , t 2 )Q(t 2 - h) 

D--{n, h; f 2 , t 2 ) = D>(fi, *i; f 2 , t 2 )<d(t 2 - h) + D<(f u * l5 f 2 , t 2 )6(ti - * 2 ) 

D+-(fi,*i;f 2 ,* 2 ) = -£><(fi,ti;f 2 ,t 2 ) 

£>- + (ri,t a ;f 2 ,t 2 ) = -D>(f 1 ,t 1 ;f 2 ,t 2 )=- J D<(f 1 ,t 1 ;f 2 ,t 2 ) (47) 

As usual, the path integral over the quadratic form may be evaluated and we may then 
expand all interaction terms perturbatively for weak couplings. 

After some straightforward algebra we obtain the following partition function 

zvt, Ji, 4] = ex P {4 dt[H jl r _ exp {4 dtH j^} 

?[-iTl+> i TT=\ xex p|-^ / dt ^ i dt 2 J%(jt 1 )D$(t 1 ,t2)Jg(jt2)\ 
0J& dJQ 1 2 Jc Jc ) 

(48) 

where J2 stands for the source term on the contour C. 
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In the limit T' — > oo and, T — > — oo the contributions from the terms in which one of 
the currents is j£ or and the other is a J$ vanish when computing correlation functions 
in which the external legs are at finite real time, as a consequence of the Riemann-Lebesgue 
lemma p0|. For this real time correlation functions, there is no contribution from the J$ 
terms that cancel between numerator and denominator. Then the calculation of finite real 
time correlation functions, the generating functional simplifies to 



5 5 



8Jf SJq, 



--/ dtj dt 2 J°(ti)Dl b (t u t 2 )JUt2)\ (49) 

Z J— oo J —oo ) 

with a, b = +, — . 

The Green's functions are similar to (46j47|) and obtained from them by replacing 
m a ; f3 a by the values for the <3> field. 

Appendix B: Evolution equations 

The most convenient method to obtain the evolution equations is the tadpole method [50|. 



In the present non-equilibrium situation, it is implemented in the following manner. First 
we recall that the non-equilibrium generating function without sources is given by 

Z[0] = J £)$i£)$ 2 £)$3 / V$ + V$-V^ e iS eff [$+,*-, $0] ( 50 ) 

with the boundary conditions $+(T) = $^(T - ifc) = $i;$~(T) = $^(T) = $ 2 and 
$+(T") = $~(X") = $ 3 . After integrating out the degrees of freedom of field a (bath), «S e // 
including the influence functional is thus 



s eff [<s> + ,<s>-,^] = J t dt{c^ + ]-c^-}} + j T l *dt 

\g 2 [ T dh [ T dt^^Df^MWth) (51) 
I Jo Jo 



We shift the $ ± fields by 

$ ± (f,t) = <f>(t) + <S ± (r,t) (52) 



where <p is a background mean field. The tadpole method |3(J requires that 

($ ± (f,t)) = (53) 
We also need to shift the field $^ by 

t) = 0(T) + ¥{r, t) (54) 
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Prior to the time when the system is coupled to the bath, the system is in equilibrium and 
is a constant in time. So for t < , <p(t) = 0(0) which is the initial equilibrium value of 
the mean field. 
Now we expand 

£ $ = C + — $ + -— $ 2 + TTTTyjvT + • • • 

5$ 2 5<3> 2 3!5$ 3 

for the (+,—,/?) branches and consider the linear, cubic and quartic terms in $ as 
perturbations. The sources are now coupled to the fluctuation part of the field to generate 
the perturbative expansion and require that order by order in perturbation theory the tadpole 



condition (|53"D is fulfilled. 



An important ingredient is the identity 

D ++ + d + - + D~+ + D = 
for both the a and $ Green's functions. 
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Figure Captions: 

Figure 1: Contour in the complex time plane for the non-equilibrium generating func- 
tional. 

Figure 2: a) Propagators (++) for the "system" field (straight line) and "bath" field 

(wavy line), there are four of each (++ ; H — ; — h ; ). b) One loop correction to 0(\g 2 ). 

There are four terms for the insertion of the bath propagators. 

Figure 3: Numerical evolution of the effective equation of motion for x(s) vs. s for the 
values Hr = m a ; g = ml ; T a = T c = 10m a . The initial conditions are: x(0) = 1 ; x(0) = 0. 
The solid line is the evolution of the classical field equations (without the non-Markovian 
kernel). The long dashed line is the evolution including the non-Markovian kernel but 
without the one-loop correction. The short dashed line is the full equation including the 
non-Markovian kernel and the one-loop contribution. 
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